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HOLOGRAPHIC MEASUREMENT OF WAVE 
PROPAGATION IN AXI- SYMMETRIC SHELLS 


David A. Evens en 
Robert Apr ah ami an 
Jerold L. Jacoby 

TRW Systems Group 


SUMMARY 

Holography is a lens less imaging technique which allows the 
reconstruction of three dimensional images. A related technique, 
called holographic interferometry, allows one to measure static or 
dynamic displacement on the order of a wavelength of light. The 
present report deals with the use of pulsed, double-exposure holographic 
interferometry to record the propagation of transverse waves in thin- 
walled axi-symmetric shells. 

The report is subdivided into sections dealing with 

(i) Wave propagation in circular cylindrical shells 

(ii) Wave propagation past cut-outs and stiffeners 

and 

(iii) Wave propagation in conical shells 

Several interf erograms are presented herein which show the waves 

reflecting from the shell boundaries, from cut-outs, and from stiffening 

rings. The initial response of the shell was nearly axi-symmetric 

in all cases, but non- symmetric modes soon appeared in the radial response 

(i.e., w w + w cos n 6). This result suggests that the axi-symmetric 
o n 

response of the shell may be dynamically unstable, and thus may 
preferentially excite certain circumferential harmonics through 
parametric excitation. 

Attempts were made throughout this study to correlate the 
experimental data with analysis. For the most part, good agreement 
between theory and experiment was obtained. Occasional differences 
were attributed primarily to simplifying assumptions used in the analysis. 
As faster and more economical computer analyses of wave propagation 
become available, it is hoped that the experimental data given herein 
will serve to verify computer solutions. 

From the standpoint of engineering applications, it is clear that 
pulsed laser holography can be used to obtain quantitative engineering 
data. Areas of dynamic stress concentration, stress concentration 
factors, local anomalies, etc., can be readily determined by holography. 
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1.0 INTRODUCTION 


Wave-propagation in thin walled shells is of interest in many 
areas of applied mechanics. For example, pyrotechnic shock loading 
in launch vehicles causes waves to propagate along the length of the 
structure. Other practical areas of application include waves interacting 
with cut-outs and stiffeners in re-entry vehicles, spacecraft, etc- 


From an analytical standpoint, the problem of calculating the 
shell response is very difficult (unless simplifying assumptions are 
made). For example, one common approximation is to assume the response 
is axi-symmetric. This assumption can then be verified by experiment, 
if necessary. Another common assumption is to neglect rotary inertia 
and transverse shear effects in the analysis. Again, this approximation 
can be checked (by using an improved theory) or by experiments. 

The primary objective of the present study is to provide 
experimental data for comparison with analysis. Since the experiments 
involve the use of holography (Refs. 1-8), the reader who is unfamiliar 
with the method may wish to look up either the texts or reports 
referenced herein. It is safe to say that by now holography is a 
widely-accepted measuring technique. Pulsed laser holography has been 
used in previous studies of wave propagation (Refs. 7, 8) in 

structures, as well as measuring shock waves in gases and solids. 

Regarding wave propagation in shells, it is worth noting that 
most of the analyses are (1) linear and (2) restricted to shells of 
revolution. A further simplification that is ordinarily introduced 
in the calculations is to ignore all variations in 0, the circumferential 
co-ordinate. Only recently have non-axis 3 nnmetric computer analyses 
become available. Thus, one objective of this study was to provide 
non-symmetric data (showing waves passing a cut-out, etc.) to compare 
with computer simulations. 

The tests were done on two aluminum cylinders, with a nominal 
radius a = 5", and the thickness of h = .106, and h = .25 inches, 
respectively. The experimental results suggest, that the axi-symmetric 
response may not be stable . That is, if you have an axi-symmetric 
wave propagating longitudinally along the shell, will it then cause 
the non-symmetric modes (cos n 6) to become unstable. The experimental 
results suggest that the axi-symmetric response is indeed unstable, 
since certain harmonics (cos 120, for example) are predominant in the 
non-symmetric response. 

Another point of interest involves the response of the thick-walled 
shell (h = .25 inches). In this case, the loading of T = 10 ysec 
duration excited an elasticity type response which could not be 
predicted using shell theory. 
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Results given herein include the holograms (photos) and in many 
cases the corresponding radial displacement, w. For the most part, 
the theory and experiment are in agreement, with the lack of agreement 
thought to be attributed to simplifying assumptions of the analysis. 

Tests were also conducted on two conical shell specimens as well, 
and these results are included herein. 

From the standpoint of practical engineering applications, it is 
clear that pulsed laser holography can be used to obtain quantitative 
engineering data. Areas of dynamic stress concentration, stress 
concentration factors, local anomalies, etc., can be readily determined 
and evaluated by holography. 



2.0 DOUBLE-EXPOSURE HOLOGRAPHIC INTERFEROMETRY 


Holographic interferometry is an optical technique which allows measure- 
ment of displacements on the order of a wavelength of light. Textbooks are 
available which discxiss holography and its applications (Refs. 1, 2, and 3). 

A brief summary of holography which is oriented toward structures and dynamics 
applications is given in Ref. 4. Detailed descriptions of holographic inter- 
ferometry are thus readily available to the \ininitiated reader. 


The present study deals exclusively with the application of double- 
exposure holography, using a pulsed ruby laser. First a hologram is made of 
the undeformed, stationairy object (in this case a shell). Then the object 
is deformed (by an applied load) and the hologram is exposed for a second 
time. When this "double-exposed hologram" is developed and then illuminated, 
two images are produced: one is from the undeformed object and the other 

from the deformed object. The light waves (which form the two images) 
interact with one another and create visible interference fringe patterns. 

By analyzing these fringe patterns, one can determine the surface deforma- 
tions of the object (see Ref. 5), 


Each exposure of the hologram by the pulsed ruby laser involves a 
light pulse which has a duration of about 50 nanoseconds (50 x 10“ sec.). 
This exposure time is sufficiently short to effectively "freeze" the deforma- 
tion of the structure at a particular instant in time. Thus the double- 
exposure hologram gives a record of the deformation on the surface of the 
object, at a particular time , say t^. This spatial recording can be used to 
complement conventional strain gages or accelerometers which yield time- 
histories of the deformation, at a fixed location , say "x^. 


Double-exposure holographic interferometry has been used to measure 
wave propagation in rods (Ref, 6), beams (Ref. 7), plates (Ref. 8), as 
well as surface waves, and shock waves in solids. Aerodynamic waves have 
also been studied using holography (Refs. 9, 10). Thus it is clear that 
holographic interferometry is now a well-established experimental technique 
which has been successfully applied to many problems in engineering 
mechanics . 

Typical double-exposure holographic results are given in Figure 1, 
which shows wave propagation in a flat plate. 
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(c) A = 130 usec 


(d) A = 150 usec 


Figure 1(a): Flexural waves in an elastic plate, recorded holographically 

(see Ref. 8) . 
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3.0 TRANSVERSE WAVES IN CIRCULAR CYLINDRICAL SHELLS 


Introduction and Description of the Problem 

Consider a thin-walled circular cylindrical shell, which is subjected 
to an axi-symmetric impulse loading in the radial direction (see Fig. 2). 




Figure 2: Radial pulse loading on shell 

(r, 6, X, cylindrical polar co-ordinate) 


The shell is described in terms of cylindrical co-ordinates (r, 0, x) in 
the radial, circumferential and longitudinal directions, respectively. 
The corresponding elastic displacements of the shell are w (radial) , 

V (circumferential) , and u (longitudinal) . The axi-symmetric impulsive 
load causes waves to propagate longitudinally, and the experiments 
described herein give results for the radial displacement (w) as a 
function of position (x) along the shell. 

The following sections discuss some analytical considerations 
relating to this problem, and experimental holographic results involving 
both symmetric and non-symmetric responses, reflections from the ends 
of the shell, and comparisons between theory and experiment. 
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Some Analytical Considerations 

For dynamic problems (such as wave propagation) in thin-walled shells, 
the mathematical modeling of the structure usually involves one of the 
following levels of sophistication: 

(I) - Lovers First Approximation (Ref. 11) which is a shell theory 
analogous to Bemoulli-Euler Beam Theory. 

(II) - "Timoshenko" Theory, which is a shell theory that includes 
rotary inertia and transverse shear effects (See Ref. 12) 

(III) - Three-Dimensional Elasticity Theory (Ref. 13), which the 
shell theories approximate. 

The main differences between these three mathematical approaches to 
shell dynamics can be described in terms of the vibration frequencies and 
structural wavelengths which are excited by the applied loading. 

For example, reference 14 indicates that conventional, thin shell 
theory (I) applies when 


0 < 0)- < .1 03 
~ f — s 


(3-1) 


where designates 


the flexural vibration frequencies and 



2(l+v) 


1/2 


(3-2) 


is the frequency of the first "thickness-shear" mode in an infinite flat 
plate (Ref. 15). 

When rotary inertia and shear are included in the analysis, the corre- 
sponding inequality is (Ref. 14) 


0 < 03^.^ < 03 (3-3) 

— II — s 

where 03 ^^ denotes the vibration frequencies computed using this ’Improved" 
shell theory. 

Elasticity theory (III) is valid for 


0 < 


03 < 0)e 


(3-4) 
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where o) is a cut-’off frequency significantly greater than o) . For example, 
may oe associated with lattice vibrations or other effects in which the 
equations of continuum mechanics no longer apply. 

In the experiments reported herein, the forcing function applied to the 
shell has a time-history given approximately by 


f(t) = sin^ ^ 


0 < t < T 


which is shown in Figure 3. The corresponding Fourier transform. 


oo 

F(,ui) = J' e sin^ ^ dt 

— OO 


(tabulated in Ref. 16) is readily evaluated and has a spectrum |f(w)| as 
indicated in Figure 4. This spectrum (of the input forcing function) has 
a well-defined cutoff frequency given by 


CO 

c 


47T 

T 


(3-5) 


where T is the duration of the input pulse. 

The experimental loading device (described in Appendix A) gave a pulse 
duration on the order of 10 psec. From Equation (3-5), the corresponding 
cutoff frequency is 

CO = 4 7 T X 10^ = 1.256 X 10^ rad/sec (3-6) 

c 


That is, the input forcing function contains frequencies primarily within 
the range 


0 < CO < 


1.256 X 10® 


rad/sec 


For comparison purposes, the thickness-shear frequency of the thin- 
walled cylinder used in our experiments is given by 


0 ) = = 3.83 X 10^ rad/sec (3-7) 

s n 

(i.e., an aluminum shell with a thickness h = .1 inch). Equations (3-6) and 
(3-7) show that in our thin-shell experiments, we have frequencies in the 
range 

0 < CO < . 3 CO (3-8) 

— — s 
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Radial pressure loading, p/p 



Non-dimensional time, (t/T) 


Figure 3: Time-history of the applied load (schematic) 
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(a) Pressure-time history 



(b) Spectrum of the pulse 


Figure 4: Input pulse and its Fourier transform 
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which suggests that conventional shell theory (I) may not apply, in view of 
the inequality (3-1) • 

A similar calculation for the thick-walled aluminum test cylinder 
(he .25 inches) gives 


w = 1.53 X 10^ 
s 


rad/sec 


and 


0 < 0 ) < ,8 0 ) 
~ — s 


(3-9) 


which clearly exceeds the bounds of the inequality (3-1) given by Ref. 14 
for conventional shell theory (I). 

However, it appears that the frequency bounds (3-1) and (3-3) given by 
Reference 14 are somewhat conservative. For example, it is possible to 
exceed these bounds (for transient loadings) and still obtain adequate 
solutions for the radial displacement, w(x,t). In such cases, conventional 
theory (I) may correctly predict the displacements but not the accelerations 
(and/or strains), which involve derivatives of the displacement. A thorough 
discussion of this point and the limitations of shell theory (I and II) is 
given in Ref. 17. 

It should be noted that the thickness- shear frequency (w ) given in 
Equation (3-2) is directly related to the time required for a shear wave to 
propagate through the thickness of the shell. That is. 


T 

s 



(3-10) 


is the travel time for a wave (having speed c ) to propagate a distance 2h 
(round trip). For the cylindrical shells which were tested. Equation (3-10) 
gives 

= 1.65 ysec (for h = .1 inch) 

and T =4.1 ysec (for h = . 25 inch) 

s 


These travel times are significantly less than the duration of the input 
force pulse (T « 10 ysec) and indicate that wave propagation through the 
thickness of the shell is much less important than wave propagation along 
the length of the shell. This gives additional evidence to support the con- 
clusion that shell theory will adequately describe the observed experimental 
response, and that it is not necessary to use elasticity theory (III) in the 
analysis. 
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With this preliminary discussion as a backgrotind, the next section 
describes the tests on a thin-walled aluminum shell and the corresponding 
holographic measurements. 
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Experimental Set-Up and Test Procedure 


Test Specimens - The shell specimens used in these experiments 
were sections of seamless aluminum tubing, which is available commercially. 
The tubing was cut to the nominal lencth of 24 inches and was of approxi- 
mately 10 inches in diameter. Two wall thicknesses were used, namely 


h = .1 inch (referred to as the "thin” shell) 


and 


h ^ .25 inch (referred to as the "thick" shell) 

The wall thickness and the inside diameter of the thin shell were 
carefully measured at several locations around the circumference. These 
measurements are summarized in Table I. 


TABLE I 


Geometry of Thin-Walled Cylinder 


Average wall thickness. 

Maximum deviation in wall 
thickness. 

Per cent deviation in wall 
thickness. 


h = 0.1060 inches 
Ah = 0.0053 inches 



Average inside diameter, = 9.807 inches 

Maximum variation in 

diameter AD = 0.080 inches 


Per cent variation in 
diameter 

Average radius to the mid- 
surface 

Radius /thickness ratio 


D. 

1 

R 

R/h 


0.82% 



46.8 


= 4.957 in 


The shell specimen was supported by a heavy steel ring which was 
pressed into the left end of the shell. (A photograph of the thin 
shell specimen is given in Figure 5) • The steel end rings were used 

to simulate a clamped boundary condition (w = 0, ^ =0, u = 0) at 

the ends of the shell. For some tests, the right-hand ring was 
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removed to provide a stress-free boundary. The impulsive load was 
applied midway along the length of the shell (see Figure 2). A 
discussion of the applied load is given in the section which follows. 

The Applied Loading - The load applied to the shell was subject 
to the following design constraints: 

(1) The load had to be repeatable , so that the same forcing 
function could be applied several times. 

(2) The timing had to be precise (i.e., to within a micro-second) 
so that the laser firing could be properly controlled. 

(3) The load had to be concentrated in space and of short duration 
in time. 

(4) The load should excite axi-symmetric waves, for comparison with 
available analyses. 

The method selected to meet these criteria involved using a 
repelling wire technique, which is described in detail in Appendix A. 

The technique relies on the fact that two parallel wires which 
carry current exert an electromagnetic force on each other. By sending 
a large surge of current through two concentric loops of wire, the 
outer wire was driven radially outward to exert a force pulse on the 
inside of the cylindrical shell. 

The force on the shell can be described (approximately) by a 
band of pressure that is 

o uniform in the axial direction (x) with a length of ^ inch 

o nearly uniform around the circumference (0), 

o of short duration in time, on the order of 10 ysec 
Mathematically, the band of pressure p(x,0,t) is approximated by 

p(x,e,t) = A sin — 


within the domain 

'E ^ X < e 

0 1 e ^ 2 tt 

0 1 t ^ T 
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and zero elsewhere. The value of the constant e is 1/8 of an inch, 
and the load duration T is approximately 10 ps. A more detailed 
discussion of the applied loading is given in Appendix A. 

Optical Arrangement for Recording Holograms - The optical arrange-- 
ment (shown schematically in Figure 6) was used to record the double- 
exposure holograms. The actual equipment used in the experiments is 
shown in Figure 7. The equipment and optical set-up are standard 
for pulsed-laser holography (e.g. , see Ref. 7). 

Timing of the Laser Pulse(s) - The holograms were made by first 
firing the laser to record the static, underformed position of the 
shell and then firing the laser a second time to record the wave 
propagation event- The initial, first pulse of the laser (first exposure) 
D resented no problems and was readily accomplished. The second laser 
pulse (second exposure) had to be synchronized with respect to the 
loading of the shell. The timing was accomplished electronically, using 
well-established techniques (Refs. 7, 8). 

An inductive loop near the capacitor bank provided the timing 
signal to mark the beginning of the repelling-wire loading pulse. This 
timing signal was used as the input to the adjustable time-delay circuit 
of an oscilloscope. The output signal from the oscilloscope was used to 
control the laser pulse, by means of the Pockels cell within the laser 
cavity. This arrangement allowed the experimenter to initiate the waves 
in the shell (by firing the repelling-wire circuit) and then to record 
the deformation at a desired time later (i.e.,. at time t = 5 ys, 10 ys, 

50 ys, etc.) by means of the laser pulse. The separation in time was 
obtained by adjusting a potentiometer on the oscilloscope. Since the 
pot settings provide only a nominal estimate of the time-delay, the 
timing signal (from the repelling-wire circuit) and the light pulse 
(which was sensed by a fast photodiode) were recorded simultaneously 
on an oscilloscope trace. This trace provided an accurate measurement 
of the timing. 

Test Procedure - Once the experiment was arranged and the electronic 
timing established, the test procedure was straightforward, as follows: 

(1) The hologram was placed in position and the laser fired to 
give the first exposure. 

(2) The capacitors for the repelling-wire circuit were charged 
to the desired voltage. 

(3) The oscilloscope time-delay was set to the desired value 
(e.g. , 50 ys) 

(4) The repelling-wire circuit was fired, which initiated the 
timing sequence and resulted in the second exposure of the 
hologram. 
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Hologram 

Plate 





Figure 7 : Experimental Set-up and Optical Arrangement 
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Typical holograms and results obtained in this fashion are 
discussed in the pages which follow* 


Typical Results for the Thin-Walled Cylinder 

Demonstration of Repeatability - Holograms were recorded at several 
different time-delays, which showed the waves at various distances along 
the shell. To record the waves at times t “ 5ys, 10 ys, etc. required a 
separate experiment for each time. That is, the shell was loaded once and 
a hologram recorded at t = 5 ys; then the shell was loaded again and 
another hologram made at t = 10 ys, and so forth. 

For purposes of understanding the wave propagation phenomena, it is 
desirable to relate the results obtained from several experiments (at 
individual times t*, where i = 1,2,3,...N) to a single experiment (in which 
data is obtained a^ various times t^) . 

In particular, the assumption is made that the holograms obtained 
from several separate experiments are equivalent to making several holo- 
grams in rapid succession on a single experiment. This assumption is 
akin to the "ergodic hypothesis” in random process theory (See Ref. 18, 
for example) in which statistics obtained from a single record (in time) 
are related to several similar records (the ensemble) . 

If the experiment is completely deterministic, and contains no 
random behavior, then holograms obtained from separate tests but at the 
same time-delay (say t ~ 50 ys) will show identical fringe patterns. Such 
behavior is indicated in Figure 8, which shows the results of three 
tests, all run for the same loading input and the same time-delay. The 
fact that these fringe patterns (Figure 8) are virtually identical 
demonstrates the repeatability of the experiment and serves to justify the 
equivalence assumption just discussed. 

Determination of Radial Displacements from the Inter ferograms - 
Figures 9 through 11 shovr interferograms obtained at the times t = 5 ys, 
10 ys, and 20 ys after initiation of the impulse. The corresponding 
curves of radial displacement (w) versus x (the distance along the shell 
length, measured from the plane of loading) are given in Figure 12. In 
order to understand how the data points in Figure 12 were obtained from 
the photographs, it is first necessary to discuss the equations relating 
the shell displacements and the interference fringes. 

The general relationship between a displacement field and the corre- 
sponding holographic fringe pattern has been examined by several authors 
(e.g. , Refs. 5 and 19 ) • In particular. Reference 19 shows that the 
reconstructed image of a hologram exposed twice to a body which has 
deformed between exposures will exhibit dark fringes on parts of the body 
wherever the condition 
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Figure 8(a): Interferogram Showing Waves at 5 = 50*5 ysec 

(Shot Number 2) 



Figure 8(b) : Interferogram Showing Waves at t = 50 ysec 

(Shot Number 3) 



Figure 8(c) : Interferogram Shox-7ing Waves at t - 50 ysec 

(Shot Number 8) 
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Figure 12: Radial Deflection Along the Shell 








(2nirl)X 

2 


(3-11) 


is satisfied; where 

the displacement vector 

wavelength of the light from the laser used to 
make and reconstruct the holographic images 

unit vector in the direction from the object 
to the illuminating source 

unit vector in the direction of view, from the 
object through the hologram to the observer 

integer, the fringe order, jW, +2, +3, etc. 

n ) is illustrated by the vector diagram shox-m in 
Figure 13 • The term (n^ + n^) is a vector sum, and it is represented 
by a vector xvrhich lies in the plane containing and n^ and bisects 

the angle betx^een them. The magnitude of this vector sum is 2 cos i(n , n ), 

->■ -> ^ i V 

X7here (u. , n ) is the angle which is bisected. Finally, the dot product 

d % (n. + n ) is the projection of the displacement vector 6 in the 

direction of (n^ + shox-m in Figure 13. 

A particular case of Equation (3-11) x^hich often occurs in transverse 
x^ave propagation problems is illustrated in Figure 14. Referring to 
Figure 14, the structure of interest (in our case, the circular shell) lies 
Xirith its axis in the horizontal plane. The angle 0 is measured positive 

clockx^ise from the horizontal plane as shoxm in Figure 14 . A general 
“>■ 

displacement, 6 , of the shell surface is represented by 


6 = 
A = 



n = 

The term 6- (n. + 


6 = 


-> 
u e 

X 


+ v 




+ X7 


e 


x^rhere e^, e^ , and e^ represent unit base vectors, and u, v, X 7 are the 

corresponding displacements, in the (x, 0, r) directions, respectively. 

The illumination and viewing vectors are given (in Cartesian co-ordinates) 
by 



+ n 


iy 


->■ 


e 

y 





n 


vx 


e + 

X 
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Observer 
or Camera 



Object 

(e.g., Cylindrical shell) 


N3 

U1 


Figure 13: Vector diagram shov/ing the illuminating and viewing directions, 

->■ 

with the displacement 6 between exposures. 



Figure 14: Horizontal plane, radial displacement, and unit vector (s) for analyzing the 

interf erograms . 


where e is the same unit vector (along the axis of the shell) and 
e is a unit vector in the (horizontal) x-y plane. The relationship 

“V- 

between e and e , e^ is 
y r * 0 


e = e cos 0 - e. sin 0 
y r 0 

-y- 

Forming the dot product 6 • + "S^) gives 


6*(n.+n)=u(n. +n ) 

1 V XX vx 


- v sin 0(n. + n ) 

xy vy 


+ w cos 0 (n . + n ) 

ly vy^ 


( 3 - 12 ) 


At this point, it is convenient to introduce two engineering 
approximations : 

(i) For transverse wave propagation, involving primarily 
flexural motions, calculations (as well as previous 
experiments) show that the radial deformation is 
predominant : 

That is. 


u<< w 
v<^ w 


(ii) In setting up the holographic experiment, it is often 
possible to locate the_^illumination and viewing points 
such that the vectors n^ and n^ are nearly perpendicular 

to the x-axis ; That is 


n. < < n. 
XX xy 


« 


vx 


n 

vy 


-y -y -y 

the vectors and are nearly parallel to e^. 
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Under these conditions 


Equation (3-12) is well-approximated by 


6 • (n^ + n^) 


w cos 0(n. + n ) 

ly vy 


or 


(n. + n ) 

X V 


= 2 cos 6 


Since the components n^^ and n^^ are approximately unity. Thus, we 
have (from Equations 3-11 and 3-12) 


(n. 


+ n ) = j 

V 2 


^ 2 w cos 0 


or 


^ (2n ± 1)X 
4 cos 0 


(3-13) 


The most direct approach to reducing the holographic fringe data is to 
apply Equation (3-13) along the horizontal plane (i.e., centerline) 
through the shell (where 0=0). Then Equation (3-13) gives 


^ (2n ± 1)X 


(3-14) 


for the radial deflection, w. 

In order to apply the preceding equations and obtain quantitative 
data, it is necessary to determine the fringe order, n, corresponding to 
each fringe seen in the reconstructed interf erograms (Figures 8 through 
11). The assignment of fringe order requires some engineering judgment 
and was influenced by the following considerations: 

(a) In response to the applied loading (which acts radially outward) the 
portion of the shell directly over the load will initially move outward, 
in the radial direction. It is expected that the maximum radial deflec- 
tion will occur at the plane of loading at least for the first few 
holograms (which record the initial response) . 

(b) Ahead of the wavefront, where no fringes appear, the deflection of 
the shell is known to be "zero”, which provides a reference level. (More 
precisely, the absence of a fringe in the "undisturbed" region ahead of 
the wave indicates that the deflection there is less than one-quarter of 
a wavelength of the laser light.) 
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(c) From previous experience with wave propagation phenomena (e.g., waves 
in beams, plates, etc.)> it is expected that the deflections will describe 
a series of ripples or undulations which are alternately up and down. That 
is, the shell surface will be a series of hills and valleys (as a function 
of the longitudinal coordinate, x) much like a damped sine wave. 

As an example of the data reduction process, consider Figure 10, 

(t = 10 ys). Starting at the center of the fringe pattern (i.e., at x = 0) 
and proceeding toward the right half of the photo, the interference 
fringes indicate a central "hill" and then a slope downward to the right. 
(The interference fringes may be regarded as contour lines on a topographic 
map. The closely-spaced fringes for x > 0 indicate that the deflection 
surface slopes downward to the right in that region.) Continuing to the 
right, the photograph shows a wide white area followed by a heavy, dark 
fringe. Close examination of the dark fringe (on the right) shows that 
it merges with the adjacent fringe (on the left) near the top and bottom 
of the shell (i.e., near 0 = + 90°). The merging of adjacent fringes at 
the extremes of the shell is due to the "cos6 effect" described mathe- 
matically in equation (3-13). 

The fact that the two fringes coincide near 0 = 90° indicates that 
these fringes have the same fringe order, n. The wide spacing between 
the two fringes at 0 = 0 suggests that they represent opposite sides of 
a "valley" in the deflection curve. To the right of this valley, no 
additional fringes are evident, which indicates a lack of deformation in 
this region. If we now retrace our steps and move from right to left 
across the fringe pattern, we conclude that the heavy dark fringe (on 
the right) has a fringe order n = -1, which corresponds to a radial dis- 
placement w = - A/4 (inward) at that poi.nt on the shell. Proceeding to 
the left (across the valley) we come to the "companion" fringe, which 
also has order n = -1 and deflection w = -A/4 (inward). Continuing to 
the left, the fringe order increases monotonically (n = 1, 2 , 3, . - . etc. ) 
as we climb the hill in the deflection surface up to the crest. 

Having thus identified the fringe order, the radial displacement 
can be computed at each fringe location using equation (3-14) . In order 
to plot the displacement (as a function of x) the location of each fringe 
was determined by mechanical measurement on the photograph using a 
traveling microscope. The curves of radial displacement vs position 
shown in Figure 12 were obtained in this fashion. 


Non-Symmetric Behavior 

Initial Response - From a qualitative standpoint, the interferograms 
shown in Figures 9 through 11 exhibit two major trends: 

(i) At the front of the wave, the interference fringes show only slight 
variation in the circumferential direction. Slight variations with 0 are 
indicated by the "cos0 effect" in equation (3-13) and indicate that the 
shell deformation is very nearly axisymmetric in the vicinity of the wave- 
front. 
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(ii) At the loading plane (i.e*, x = 0) significant circumferential varia- 
tion is indicated by the closed elliptical fringe contours. This non- 
symmetric response is evident even at the very early times (e.g., t = 5 ps) , 
including the times when the loading occurs (0 t _< 10 ps) . 

Since the applied load (due to the repelling wire) did not extend fully 
around the circumference, it was expected that non-symmetric behavior would 
eventually occur. However, the "gap" in the loading would cause waves to 
propagate circumferentially and they are required to travel half way around 
the shell to reach the front side (at 0=0). The travel time of these 
circumferential waves is given approximately by 


t 


1 


ttR _ 7t( 5 in ) 

c .2 in/ps 
P 


^ 75 psec 


The fact that non-symmetric behavior occurs much sooner than this (e.g., 
at t = 10 ps) points out that the applied loading must contain slight local 
non-uniformities (in addition to the "gap" region just discussed). 

Additional Results and a Computed Fringe Pattern - Additional results 
which exhibit significant non- S 3 mime trie behavior are the interf erograms 
shown in Figures 15 and 16. Referring to Figure 15, (which was recorded at 
time t = 75 ys) the appearance of the closed, nearly-circular fringes and 
their repeated nature (around the circumference) suggest a radial deflec- 
tion on the form 


w(x,6) = A^f^(x) + A^f^(x) cosk0 


(3-15) 


Equation (3-15) represents an axisymmetric response combined with a non- 
symmetric mode. The value of the harmonic k and the relative amplitudes 
and A^are of interest- As a rough estimate, the functional form 
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w(x,0) = [A^ + A^ cosk0] cos— ^ 
was used to approximate the response near the loading plane (x = 0) . 


(3-16) 


Using the approximations discussed previously, we have (see Equation 
(3-13) 


(2n-l)A 
4 COS0 


(3-17) 
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Figure 15 : Interf erogram Showing Waves in the Thin Shell at 

t = 75 ysec. (Note: Non-symmetric response, 

which involves cos k 6) 


where 6 is the circumferential coordinate. Combining equations (3-16) and 
^3-17) gives 


A 

o 


[1 + r coskG] CO 



(2n-l)X 
4 COS0 


(3-18) 


where r = is the ratio of amplitudes. 

Equation (3-18) was solved numerically to give the fringe contours 
X = x(6) for various values of n, the fringe order. The results are shown 
in Figure 17, where x/L is plotted vertically (for the range 0 j< x/L ^ 1/4) 
and 0 is the abscissa (0 _< 0 ^ tt/ 2). The contour lines of Figure 17 
vary from n = 1 (the outermost contour) to n = 10 (a single point at 
X = 0, 0=0). Figure 17 was calculated for k = 12, A = 12(X/4) , and 
A. = 7(X/4) and was meant to correspond with the interf erogram of Figure 15. 
Thus, the closed, circular contours of this interf erogram suggest the 
presence of cos (120) in the shell response. 

Displacement Curves - Figure 16 (recorded at t = 99 ys) shows non- 
symmetric behavior, which is also thought to involve cos(126). The shell 
response at these later times is clearly two-dimensional, involving varia- 
tions in both the longitudinal (x) and circumferential (0) directions. A 
complete description of the experimental data would involve longitudinal 
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Figure 16: Interferogram showing waves in the thin 

shell at t = 99 ysec. (Note: non-symmetric 

response, which involves cos k 0,) 
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z(SHELL AXIS) 



0 (CIRCUMFERENCE) 


Figure 17: Computer-Generated Fringe Pattern 

(See Equation 3-18) 



traverses ("cuts" at constant 0) as well as circumferential (cuts at 

constant x) . This detailed description of non-symmetric responses is a 

formidable task in data reduction and was not attempted. Representative 

plots of the radial displacement vs. x are shown in Figures 18 and 19, 

which were obtained by interpreting the fringes along the centerline 

(0=0). The experimental data shown in these figures is typical of 

flexural waves propagating along the shell in the x-direction. (See Appendix B.) 

The reader will note that the data points (see Figure 18) give a 
good definition of the shell behavior near x = 0, where several fringes 
occur within one structural wavelength of the deformation. Near the front 
of the wave, however, the wave amplitude is reduced, and only one or two 
fringes occur to characterize the deformation. Thus, for these later 
times, the oscillations near the front of the wave are not well-defined 
by the experimental data. This lack of definition at the wavefront can 
be overcome by increasing the impulse to the shell, i.e., causing an 
increase in the structural response. Increasing the shell response will 
lead to other problems, however, such as very high fringe density and 
possible lack of resolution near the plane of the loading, (x = 0) . The 
repelling wire loading device was designed! using preliminary estimates 
for the shell behavior, and consequently it was not capable of prodlucing 
large- amplitude response near the wavefront. 

Dynamic Stability of the Axi-symmetric Response - Despite attempts 
to produce a nearly axis 3 nmnetric response of the shell, the experimental 
results show that the initial axisymmetric behavior is followed by waves 
propagating in non-symmetric modes (i.e., cosk0) . The fact that one of 
these circumferential harmonics predominates suggests that the axisymmetric 
response may be dynamically unstable. In other words, an axisymmetric 
response may (through non-linear coupling terms) excite non-symmetric waves 
in the shell. In this case, small non-symmetric perturbations will initially 
increase with time, as they extract energy from the S 5 unmetric response. 

Structural behavior of this type (where a S 3 minietric response becomes 
unstable and excites cos n0 modes) is we 11- documented in maoriLy problems 
involving shells (see Ref. 20 , for example). A simple example of dynamic 
instability of this type was studied by Mclvor (Ref. 21) with reference 
to circular rings. Mclvor showed that a uniform radial imptilse excites 
the axisymmetric, "ring mode," response and that the symimetrlc response 
in turn causes certain flexural modes (involving cos nS) to participate 
in the motion. Stability studies such as Ref. 21 involve using nort-linear 
equations for the shell, whereas virtually all analyses to date for wave 
propagation in shells use linear theory (e.g. Refs. 17 , 22, and 23 ). 

From a practical standpoint, the observation of significant non- 
symmetric behavior indicates that 

(a) Calculations of axisymmetric wave propagation behavior in shells are 
rather idealized, 

(b) Considerable care must be exercised in such wave propagation experi- 
ments to eliminate circumferential variations in the structure and its 
loading, and 
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(c) Analytical studies concerning the dynamic stability of an axisymmetric 
wave might be a fruitful topic for further research. 


Waves Reflecting from the Ends of the Shell 

Background - This section deals with the propagation of transverse 
waves toward the ends of the shell, and the subsequent reflection of 
the v;ave(s) from the shell boundaries. The characteristics of the 
reflected wave depend upon the boundary conditions at the ends of the 
shell. Three different boundary conditions v/ere used, vzhich will be 
referred to as (i) "clamped’*, (ii) "free" , and (iii) "simply-supported." 
These generic terms correspond to three distinct, idealized boundary 
conditions commonly applied in theoretical studies of shell dynamics. 

The problem of an axi-symmetric wave reflecting from an axi-symmetric 
boundary is more readily analyzed than wave reflections which have a 
general, non-S 3 nnmetric character. Thus, it was desired to generate an 
axi-sjmimetric wave experimentally for comparisons with analysis. However, 
the previous test results (c.f.. Figure 15) had shovm that axi-symmetric 
behavior occurred primarily near the front of the wave, followed by 
pronounced non-S 3 mmietric response. In order to keep the reflections 
"as axi-symmetric as possible," the loading device was moved to a location 
near the boundary being studied. 

This re-adjustment of the plane of load application had another 
beneficial effect, in that the propagating wave reached the boundary 
in a shorter time. For purposes of analysis, it is more convenient 
to have the loading plane near the shell boundary, thereby reducing 
the time necessary to compute the wave reflection (i.e., shorter computer 
"run times" are required) . Despite this attempt to facilitate the 
analysis, it V7as not practical to compute the wave reflections. The 
reflection data presented in the next few paragraphs is typical of 
the experimental results for the three boundary conditions. It is 
anticipated that this data may eventually be compared with analysis when 
faster computer programs become available. 

Reflection from a Clamped Boundary - The idealized, "clamped" 
boundary (used in analysis) is specified by the conditions 

u = 0 (zero axial displacement) 
w = 0 (zero radial displacement) 

and 

^ = 0 (zero slope) 
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x^hich are applied at the end of the shell that is “clamped" . These 
conditions of zero slope and displacement were simulated experimentally 
by inserting a large steel ring in the end of the aluminum shell 
(c.f., Figure 5). 

The steel ring had a rectangular cross-section 1 inch deep and 
2 inches wide, x^hereas the aluminum shell X'7all x^ras approximately .1 
inches thick. The outer diameter of the steel ring was a fex^ thousandths 
of an inch greater than the inside diameter of the shell, and the ring 
x-jas pressed into the shell (e.g,, a force fit) to provide a x^/ell-restrained 
boundary. 

The test procedure used was similar to that described previously, but 
the loading device was located at a distance of nine (9) inches from 
the left edge of the aluminum shell. The steel end ring was inserted 
approximately one (1) inch into the end of the shell, which resulted in 
a distance L = 8 inches from the plane of the loading to the "clamped" 
boundary. 

A typical interf erogram (taken at time t = 90 ysec) is shovm in 
Figure 20. The corresponding radial deflection curve (along the line 
6 = 0) is shoxm in Figure 21. Aluminum has a shear wave speed c s:; .1 
inch/ysec, and at a time t = 90 ysec the transverse xjavefront has^ 
propagated from the loading plane, to the "clamped" boundary, and is 
starting back tov7ard the center of the shell. Several additional 
interf erograms, x^hich show the wave at various other times (before and 
after reflection) are given in Appendix C. 

Reflections from a Free Edge - A series of interferograms were made 
to determine the reflection of a transverse wave from a free edge. The 
idealized boundary conditions that the edge be "stress-free" x^ere well- 
satisfied in the experiments. Thin shell theory requires that 


and 


M = 0 

X 

(moment - free) 

o 

II 

(zero axial stress) 

II 

o 

^zero transverse shear) 


x-7hich can be related to the (more exact) boundary conditions of 
elasticity theory, namely that the surface tractions vanish: 


V 

T. 


X 


= vj .= 0 


i = 1,2,3 
j = 1,2,3 


(see Ref. 24)* The shell was cantilevered from its left end, and the 
right end x^as indeed "stress-free" in the experiment. 
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Figure 20s Wave interacting with a "clamped" boundary. (Steel ring is 
at the left of the shell. Time t = 90 psec.) 






Transverse Deflection, w, (inches) 
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j Figure 21: Radial Deflection, w, as a Function of Position along the Shell. (Wave Reflecting from a 

i fixed Boundary; see Interf erogram Figure 20.) 

i 
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A typical interferogram (showing the wave reflecting from a free 

edge) taken at time t = 105 ysec is given in Figure 22. The corresponding 

deflection data (w vs x) is plotted in Figure 23. The radial load 

was applied at x = 2.5 inches, to the right of the center of the shell. 

At time t = 105 psec, the transverse wavefront (which moves at 

approximately c ^ .1 inch/psec) has reached the edge of the shell and 
s 

Is reflecting back toward the center of the shell* 

Additional interferograms (for this same free-edge boundary condition) 
are presented in Appendix C and show the wavefront at various times. 


Reflections from a Simply-Supported Boundary - Another edge condition 
which is commonly used in analysis is the simply-supported boundary: 


= 0 (moment - free) 

w = 0 (zero radial displacement) 


The accompanying in-plane boundary condition is either 


u = 0 (zero displacement) 

or 

= 0 (zero axial force) 

An attempt was made to simulate the "simply-supported" boundary 
experimentally by machining vee grooves around the circumference of 
the shell. The grooves (i-e., "notches") produce a significant (local) 
reduction in the flexural stiffness, since 

(i) the shell thickness is reduced by a factor of 5 or 10 
at the groove 

and (ii) the bending stiffness varies as the third power of the 

3 

wall thickness, i.e., h 

Enough material is left (at the bottom of the vee) to approximate the 
restraints 

\<i ^ 0 
u « 0 

prescribed in the analysis. This is a common technique of approximating 
simply-supported edges experimentally. 
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Figure 22 : Wave interacting with a "free” edge (at the right) . Time t = 105 psec. 
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Figure 23: Radial displacement, w, as a function of axial position, x. Load applied at x » 2.5 inches, 

free-edge at x « 8.5 inches. (See interferogram, Figure 22.) 
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The vee groove was situated at the left end o£ the shell, close 
to the steel end ring. A typical interterogram (showing a wave reflecting 
from the vee groove) at time t « 39 ys is given in Figure 24. The 
corresponding deflection data (w vs x) is plotted in Figure 25. 

The impulsive load was applied at a location 4 inches to the right 
of the vee groove, as illustrated in Figure 25. Additional interferograms 
(for this boundary condition) are included in Appendix C. 


Comparison with Computer Analysis 

One objective of the experimental work was to provide data for 
comparison with computer analyses. The tests on the "thin” shell (h = .106 
inch wall thickness) were co-ordinated with NASA/Langley for comparison 
with an axi-symmetric, finite-difference analysis (Reference 25). 
Quantitative details of this comparison are presented herein. 

Input Data for the Analysis - The computer program used the following 
input data: 


Radius, a = 5.0 inches 
Thickness, h = .106 inches 

Length, L = 12.0 inches (total length between ends) 
Radius/Thickness Ratio, a/h = 47.1 
Length/Radius Ratio, L/a = 2.4 

The applied loading was modeled as a pressure, p(x, t) as follows: 

[ 2 3 4 1 

16 (|) - 32 (|) + 16(|) Jh (3-19) 

for 0 ^ t ^ T (the pulse duration) 

in a region 1/4 - of an inch wide, centered at x = 6 inches. (That is, 
the load was a uniform band of pressure, located mid-way between the 
ends of the 12 inch long shell.) The factor, p, in Equation (3-19) 
represents the mass density of the shell and the constant k was used 
as a scaling factor to adjust the amplitude of the impulsive loading. 

Equation (3-19) is a polynomial representation of the applied 
pulse time-history, and it is quite similar to the expression 

p(x, 0, t) -A sin^ — 
discussed previously (c.f. Figure 3). 
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The pulse duration, T, used in the analysis was T = 10 psec, which 
closely represents the experimental result; (see Appendix A). The 
material properties used were nominal values for aluminum, namely 


E = 10 X 10^ Ib/in^ (Young Modulus) 
V = 0.3 (Poisson's Ratio) 

3 

pg = 0.10 Ib/in (Density by weight) 


The applied impulsive load, p(x,0,t) is characterized by its "specific 
impulse," which was calculated as follows: 


I 

o 



dx dt 


(3-20) 


V/hen Equation (3-19) is substituted for p(x,t) in Equation (3-20) and the 
integral is evaluated, the result is 


= ph(-i)k(|j) X lO"-^ (3-21) 

where the (1/4) comes from integrating on x from e = -1/8 to e = +1/8 
(i.e., the width of the applied load). The product (ph) represents the 
shell (mass/unit area), and the factor 8/15 x lO'-^ results from integrating 
the polynomial from t = 0 to T - 10 sec. The constant k was chosen to be 
3.65 X 10^, which yields a value 


I = 1.34 X 10 
o 


lb- sec/ in 


(3-22) 


when used in Equation (3-21) along with the nominal values of h and p. 
Thus, the scaling constant k was chosen such that the analysis used the 
same "impulse" as given by the pendulum measurements. (See Equation A-12, 
Appendix A.) 
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Ratio Between Experimental /Theoretical Responses 


With the scaling constant k adjusted in this fashion, the ratio 
bet^-7een experimental and computed deflections was expected to be 1:1* 
However, the computer solutions were noticeably lower in magnitude than 
the response measured experimentally. This result V7as disturbing, 
expecially in view of the fact that previous holographic measurements 
of wave propagation (in beams. Ref. 7, and plates. Ref. 8) had shoxm 
excellent correlation mth analysis. 

At this point, it was suspected that the pendulum calibration of the 
specific impulse might be in error (see Appendix a;. The calculation of 
the specific impulse (using electromagnetic theory and the measured 
current, I(t), in the repelling-wire device) resulted in 

2.06 X lO"^ < (I ) ^ < 2.39 X lO"^ (3-23) 

— o calc — 


where the units are Ib-sec/inch and the inequality (3-23) depends on the 
limits of an integration, discussed in Appendix A. 

T^en the value 


(I ) , = 2.39 X 10 

o calc 


-4 


lb-sec/ in 


is used to ”scale up" the computer results, (by the ratio K = Y34 

fairly good agreement is achieved between theory and experiment, as 
illustrated in Figure 26. 


1 . 8 ) 


Figure 26 shows the computed results (Ref. 25) and the holographic 
data, for times t = 5, 10, and 20 psec. As shown in Figure 26, qualita- 
tive agreement between theory and experiment is quite good, and the data 
agree quantitatively as well, providing we scale the theory to correspond 
with the impulse - 2.39 x 10“4 Ib-sec/in given in Appendix A. 

Further comparison bet^-yeen theory and experiment (for the thin-walled 
shell) was hampered by two considerations: 


(i) The computer program was based on axial symmetry, whereas 
non-symmetric responses occurred (at the later times) 
experimentally. 

(ii) The computer program involved long running times (e.g. on 

the order of 30 minutes to propagate the wave for t = 20 ysec) 
and it was not practical to compute the response at the later 
times. 


Additional comparisons between theory and experiment (but for the 
thick-walled shell) are given in Section 4.0 xyhich follows. 
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4.0 WAVE PROPAGATION PAST CUT-OUTS 
AND STIFFENERS 


Background and Summary 

The problem of axi-symmetric wave propagation in a thin-walled 
shell of revolution represents a rather idealized situation from a 
practical standpoint* This is to say, that many practical problems 
(in actual engineering applications) involve shells with cut-outs, 
stiffeners, eccentricities, etc. that depart from a mathematical model 
which requires axial symmetry. Computer programs for non-S 3 nnmetric 
wave propagation in such shells are of fairly recent development 
(e.g*. Ref. 17). One of the objectives of the present experimental 
study was to provide data for comparison with such analyses. 

Experimental data are given in this section for transverse wave 
propagation past square cut-outs, circular cut-outs, and an (axi-S 3 nnmetric) 
ring stiffener. Tests run on a thick-walled shell (h = .25 inches) 
were compared with analytical results of Ref. 17. Qualitative argeement 
between theory and experiment was obtained, and the quantitative dis- 
crepancies are largely attributed to rotary inertia and transverse shear 
effects not included in the analysis. 

Other tests on the thin-Xi7alled shell (h = .106 inches) are also 
presented herein, but analytical results for this geometry were not 
available. The "thin" wall thickness (h « .1 inch) requires a very 
"fine" finite-difference mesh in the analysis, which results in "long" 
computer run times. As faster computer analyses become available, 
efforts should be made to correlate them with the present experimental 
data. 
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Tests on a Thick-Walled Cylinder 


Background - Concurrent with our experimental work, NASA/Langley 
funded an analytical study (involving a two-dimensional, finite-difference 
approach) for nonsymmetric wave propagation in shells (See Ref. 17), It 
was understood that TRW and Lockheed personnel would coordinate these 
programs and provide each other with their results. Discussions were held 
with the NASA Technical Monitor (Dr. J. P. Raney) and with T. L. Geers 
(Lockheed, Palo Alto) until we agreed upon a shell configuration which was 
acceptable for both the holographic experiments and the computer calcula- 
tions. The present section deals with the holographic results and their 
comparison with analysis. 

Shell Specimen and Test Procedure - The cylindrical shell used in 
these tests was aluitrinuiii tubing^ and it had a nominal wall thickness 
h = .25 inches. The shell was 24 inches long and had an average inside 
diameter = 9.496 inches. Thus, the radius to the mid-surface of the 
shell was a - 4.87' (inches) and the radius/thickness ratio was a/h = 19.5. 
The shell was supported at the ends by flexible aluminum discs (1/16-in. 
thick) which were designed to provide "freely-supported" boundary condi- 
tions (see Ref. 4). (A photograph of the "thick-shell" specimen is given 
in Figure 27). 

Since the "thick shell" specimen was two- and- one- half times as thick 
as the **^thin shell"' (Section 3.0), the stiffness and mass were also 2-1/2 
times larger. Elementary response calculations show that the radial dis- 
placement for our problem is inversely proportional to the shell thickness. 
Thus, in order to maintain an adequate response (and produce sufficient 
fringes on the interf erogram) it was desirable to increase the impulse 
applied to the shell and thereby compensate for its increase in thickness. 
The repelling-wire loading technique (described in Section 3.0 and Appen- 
dix A) was used, and the impulse was increased by using (i) a higher 
charging voltage, V = 25 kilovolts, and (ii) a smaller wire spacing, 
d = 0.27 inches- Using these values in equations (A-6) and (A-7) of 
Appendix A results in a specific impulse of (approximately) 

I = 6.0 X 10 ^ lb. sec. /in (4-1 ) 

o 


applied to the thick shell. 

The optical arrangement , timing of the laser pulse, and test proce- 
dure were identical to that described previously in Section 3.0. Experiments 
Were performed first on the virgin shell (without cutouts) , and then 
rectangular cut-outs were machined in the shell and the tests repeated. 
Typical results are discussed in the following paragraphs. 

Response Without Cut-Outs - Typical holographic inter ferograms of 
the thick-walled cylinder are shown in Figures 28 through 30. Figure 28 
(recorded at t = 9 ys) illustrates the axisymmetric nature of the initial 
response. Slight nonsymmetric behavior is evident in Figure 29 (t = 15 ys) 
and Figure 30 (t = 38 ys) . The corresponding curves of radial displacement 
(w) vs axial position (x) are given in Figure 31. Additional inter ferograms 
showing the response at other times are presented in Appendix D . 
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Figure 27 


Thick-walled shell specimen (aluminum 
a = 4.87 inches, h = .25 inches). 
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Figure 28 ; Wave propagation in the thick-walled shell, (h = .25 inches) 
Time t = 9 ysec. 
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Figure 29: Wave propagation in the thick-walled shell, (h = .25 inches). 

Time t = 15 ysec. 





Figure 30: Wave propagation in the thick-walled shell, (h = .25 inches). 

ui Time t = 38 ysec. 

Ui 



Figure 31 : Radial deflection, w, vs. axial position, x. Thick-walled shell, h « .25 inches, 

without cut-outs. 



The displacement curves just discussed (Figure 31) Illustrate a 
technical difficulty encountered previously: Each structural wavelength 

of the deformation is defined by very few data points (l.e*^ fringes). 

This problem could be overcome by increasing the impulse to the shell, 
but this approach was complicated by the fact that the electrical cir- 
cuitry would not permit higher voltages (as it was, we operated well 
above the ratings of several components). An alternative means of 
increasing the Impulse was to decrease the wire spacing, d, still further 
but such an approach magnifies the circumferential variations in the 
loading, which is undesirable. Thus, the decision was made to keep the 
experimental design fixed and proceed to tests using the shell with cut- 
outs. The results of those experiments are reported in the following 
paragraphs. 

Response with Cut-Outs - The same thick-walled alumimmi cylinder 
was used for these tests, except that it was machined as shown in 
Figure 32. (A photograph of the shell with cut-outs is given in Figure 33 ). 
Eight symmetrically placed cut-outs were used, since the mathematical 
model (see Ref. 17) required several planes of symmetry to minimize the 
number of degrees-of-f reedom (i.e., to cut down on the extent of the 
finite-difference mesh). 

The holographic results are shown in Figures 34 through 38, which 
graphically illustrate the waves approaching and propagating past the cut- 
outs. Reflected waves leaving the cut-out and traveling back to the center 
of the shell are also apparent in the interferograms. Since the shell 
response is clearly two-dimensional (involving both x and 0 ) , it was 
necessary to reduce the experimental data for various "cuts** along the 
shell. This was accomplished by determining the fringe order, n, and 
calculating the corresponding displacement (as outlined previously in 
Section 3.0. Cuts were taken parallel to the axis of the cylinder, and 
primarily at two angular stations, 0=0, and 0 = tt/ 12. These values of 
0 correspond to the middle of the cut-out and the top edge of the cut-out, 
respectively. 

Plots of the radial deflection vs x are shown in Figures 39 through 
47. Until the wavefront reaches the cutouts, the response of the shell 
is basically axisymmetric (cf. Figure 34), As the wave interacts with 
the cutout, the transverse response tends to be large near the edge of 
the cutout, as might be expected (cf. Figure 37). These experimental 
results are in qualitative agreement with the computer analysis, as dis- 
cussed in the next section. 
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Section A-A 


Figure 32: Sketch showing location of cut-outs in the thick shell, 

(All dimensions are given in inches.) 




Figure 33: Photograph of thick-walled shell specimen 

with S 3 nnmetrical cut-outs (c.f. Figure 32). 
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Figure 34 ; Wave interacting with cut-outs. Thick-walled shell, 
h = .25 inches. Time t = 20 psec. 




Figure 35 : Wave interacting with cut-outs. Thick-walled shell, 

h = .25 inches. Time t - 29 psec. 
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Figure 36: Wave propagation past cut-outs in the thick-walled shell 

Time t = 41 psec. 



Figure 37 : Wave propagation past cut-outs in the thick-walled shell. 

Time t = 72 ysec. 
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Figure 


38: 


Wave propagation past cut-outs in the thick-walled shell. 
Time t = 82 psec. 
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Figure 41: Radial deflection, w, vs. axial position, x, at time t * 72 usee 

shell with cut-outs 
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Figure 42 : Radial deflection, w, bs. axial position, x, at time t = 72 usee, 

shell with cut-out. 
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Comparison with the STAR Code 


As indicated previously, one objective of the experimental work 
was to provide data for comparison with computer analyses. In particular, 
the tests on the ’’thick shell” (h = .25 Inches) were designed for comparison 
with the STAR code (Ref. 17). Quantitative details of this comparison are 
presented herein. 

Input Data for the Analysis - The shell which was analyzed used 
the following input data: 

Radius /Thickness Ratio, a/h = 20 

Length/Radius Ratio, L/a =4.8 (total length between ends) 

Conditions of S 3 nnmetry were applied on the following planes: 

• The plane x = 0, which contains the radial loading, 

• The plane 0=0, which passes through the cut-out and contains 
the cylinder axis, 

• The plane 0 = tt/ 4, which passes mid-way between the cut-outs 
and also contains the cylinder axis. 

The applied loading was modeled as a transverse shear, Q , as 
follows: ^ 



Q^(O,0,t) 


5 sin^ ~ 0<t<T=.4- 

T c 

0 t > T 


Boundary conditions of ”f ree-support” , i.e.. 


w = 0 


N 

X 


= 0 


at X = + L/2 


V = 0 


were applied at the ends of the shell. 
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Ratio Between Experlmental/Theoretical Responses - To compare 
theory and experiment, the following factors were used: 

a = 5 inches (radius) 

c = .2 in/ys (wave speed) 

L/2 = 12 inches (half length of the shell) 

E = 10 X 10^ Ib/in^ (Young’s Modulus) 

V = .3 (Poisson’s Ratio) 

Thus, T = .4 I = 10 ys, the duration of the applied loading- The 

pressure loading, p(x,t), (assumed axisymraetric) can be related to the 

shear force Q as follows: 
a 

G 

p(x,t)dx = 2Q^(0,e,t) (4-2) 

“G 



where g is a small quantity representative of the width of the applied 
load. Equation (4-2) is obtained by considering equilibrium of a 
shell element in the radial direction. 

If we now integrate both sides of this equation in time, we obtain 



t)dxdt = 


I 

o 



Q^(O,0,t)dt 


(4-3) 


where I is the "specific impulse", in Ib-sec/inch. (See Section 3.0. 
o 

Equation 3-20 , for example.) Thus, Equations (4-2) and (4-3) can be 
combined to give 


I = 
o 



I 

o 


IQEa 

(1 - V^) 


(l) 


(4-4) 
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have 


Using the nominal values of E, a, v, and T (given previously) we 


_ 10(10 X loS(5) I 10 X 10~^ 1 „ 

o ■ (.91) \ 2 / " 


75 X 10 


3 l b-sec 
in 


which characterizes the (analytical) input. 

The experimental value of the specific impulse was calculated using 
Equation (A-6) and (A-7) , following the procedures outlined in Appendix A. 
Thus, from Equation (A-7) we have 


I = 
o 


final 2 / 

/ -(^) < 5 > 


.409 


.475 


(4-5) 


-2 

where for the thick shell V = 25,000 volts and d = .691 x 10 meters, 
which gives 


I = 
o 


5.86 


6.80 


X 10 


-4 lb-sec 


in 


when expressed in English units . 

The ratio of experimental to theoretical values is 


w I 

R = ^xpt ^ expt 


t/e w - I., 

theory theory 


X 10-^ 

6.80 

2.75 X 10' 


-1 


2.14 


(2. 


48 


X 10 


-7 


The results of Ref. 17 are given in terms of deflection/radius , (i.e., 
w/a) , and in order to compare theory with experiment we have 


w 

theory 

- <i> 


(a) 





j2.l4j 





■<!> 

(2.48) 

X 10 "^ X 5 

(4- 

-6) 

w 

theory 

■<!> 

(1.125 

X 10 "^) 

(4- 

-7) 

where the values of (w/a) are 

read 

from Ref. 17. 
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Equation (4-7) states that the analytical results of Ref. 17 

should be scaled by (approximately) a factor of 10 ^ to correspond 
with the experimental input conditions. 

Analytical Results and Comparison with Experiment - Typical 
experimental data obtained holographically are shown in Figures 43 through 
45, along with analytical values which have been scaled from Reference 

17. (The scaling factor used was K = 1,125 x 10 which was used to 
multiply the analytical values of w/a. (See Equation (4-7) and the 
related discussion.) The analytical and experimental curves have many 
features in common, but certain basic discrepancies are also apparent. 

In order to understand these differences, it is worthwhile to review 
some of the limitations of thin shell theory, which is used in the 
computer analysis. 

As discussed in Section 3.0, thin shell theory is an approximation 
to the more exact linear elasticity theory. In terms of the structural 
wavelengths (say A ) involved, shell theory will apply when the ratio 
(fV /h) is sufficiently large. For example, Bernoulli-Euler type shell 
theory (Ref. 11) neglects rotary inertia and shear effects, and is 
expected to be accurate when the wavelength/ thickness ratio (A /h) is 
approximately 10 (or larger). Timoshenko type shell theory (Ref. 12) 
which includes rotary inertia, etc., is valid for shorter wavelengths, 
say A /h on the order of 2 (or perhaps 1). And finally, for yet shorter 
wavelengths, more exact theories are required. 

If we keep these limitations in mind and now refer to Figure 40 j 
we observe that the analytical curve possesses rather "high-frequency" 
(i.e., short wavelength) ripples near the front of the wave. For 
example, near x = 2 inches (Fig. 40 ) the analysis exhibits a local 
structural wavelength A ~ 1 inch. Since the thickness h = .25 inches, 
we have (locally) A /h « 4, which is outside the range of validity 
/h > 10) of the Bernoulli-Euler shell theory used in the analysis 
(Ref. 17). In other words, the analytical curve (c.f. Fig. 40) is 
not expected to be valid in the high-frequency, short-wavelength region 
near the front of the wave. Accordingly, one would not expect the 
present experimental data to agree with Ref. 17 near the wavefront. 

When rotary inertia and transverse shear are included in the 
analysis, the "high-frequency ripples" (predicted by Bernoulli-Euler 
theory) are removed (Ref. 12). The short-wavelength ripples of the 
less-accurate (Bernoulli-Euler) theory are replaced by longer-wavelength 
deflections with smaller amplitudes. These effects (of shear and rotary 
inertia) are well-documented in the literature of elastic wave propagation 
(c.f. , Ref. 26) . 

Thus, the disagreement between analysis and experiment (at the 
wavefront) is attributed primarily to the shortcomings of Bernoulli- 
Euler shell theory. Of greater concern was the lack of quantitative 
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Radial Deflection, w, (inches) 



Figure 43: Comparison of theory and experiment, at time t = 20 usee 









Radial Deflection, w, (inches) 



Figure 45: Comparison of theory and experiment (theory at 50 usee, tests at t =» 47 usee) 



agreement between theory and experiment near the center of the shell 
(i.e., X = 0) where the load was applied. In this ^center" region, the 
structural wavelengths are fairly long (at least for the later times, 
such as t ^ 20 ysec) and the analysis (Ref. 17) was originally expected 
to apply. However, discussion of this problem with several colleagues* 
has lead the authors to conclude that an analysis of the subject problem 
(using Bernoulli- Euler shell theory) may in fact be inadequate even in 
the "long-wavelength** region near x = 0. Time would not permit an 
adequate study of this problem, (from an analytical standpoint) but some 
simple **demonstrative examples** are discussed in the following section 
to support the authors* conclusions. 


Some Simple Wave Propagation Problems 
in Thick and Thin-Walled Shells 

Consider the (idealized) problem of an infinitely long, circular 
cylindrical shell with thickness h and radius, a, subjected to an 
axisymraetric impulsive ring load applied at x = 0. (See, for example 
Figure 2). The shell itself has just two characteristic lengths, namely 
h, (thickness) and a, (radius) . 

If the impulsive load is idealized as 


p(x, 0, t) = 6(x) 6(t) 


(4-8) 


where 6 (x) and 6(t) are Dirac-delta functions in space and time, 
then there is no characteristic length (or characteristic time) associated 
with the loading. Finally, suppose that we consider only a single wave 
speed, c, associated with wave propagation. 

When such an idealized problem is analyzed by Bernoulli-Euler shell 
theory, wave propagation through the thickness, h, is ignored. Thus, 
there exists only one characteristic length (namely a, the radius) in 
the Bernoulli-Euler problem, and the corresponding characteristic time 
is 


B.E c 

which is the time for a wave to travel once around the shell. In 
fact, if we use 


(4-9) 


c = c 

P 



(4-10) 


* Including Drs. T. L. Geers (Lockheed, Palo Alto), G. A. Hegemier 

(Univ. of Calif, San Diego) and J. P. Jones (Aerospace Corp., El Segundo) . 
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(see Ref. 17) for the wave speed in Equation (4-9), the result is 




where 


ring 


2u 

0) 

o 


(4-lla) 


(0 


o 


1 

a 



(4-llb) 


is the frequency of the axi-symmetric, plane strain, ”ring mode” 
(c.f.. Ref. 21). (This simplified discussion is based on the tacit 
assumption that the shell is very ”thin”, in the sense that the 
thickness/radius ratio, h/a , is very small - i.e., h/a 0.) 

Now consider the same problem when analyzed by a Timoshenko -type 
shell theory. Wave propagation through the thickness is not ignored, 
and the problem contains two characteristic lengths (h and a) with the 
associated times 


_ 2Tra _ 2 tt 

ring c c w 
p o 

(around the circumference) 
and 

^ ^ 2h ^ ^ ^ 2^ 

ts C C OJ 

s s 

where 



(4-12) 


(4-13a) 


(4-13b) 


is the lowest simple thickness- shear frequency (see Ref. 15) discussed 
in Section 3.0. (Again it is assumed that the ratio h/a is small, 
tending toward zero.) 


In writing Equations (4-12) and (4-13a) we have actually used two 


wave speeds, c and c , for the speed. 


**c" 


P s 

define two additonal characteristic times, namely 


We could equally well 


T 


s/e 


2Tra 

c 

s 


(4-14a) 
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and 


^o/h 


2h 


c 


p 


(4- 14b) 


(which together with Equations (4-12) and (4-13) make a total of four) . 
However, the latter two (Equations 4-14a and b) are of minor significance 
for the problem under discussion. 


The preceding paragraphs suggest that Bernoulli-Euler shell will 
have a time-dependence of the form 


w(x, t) ^ A cos(oj t + tf)) (4-15) 

o o 

(since is the only characteristic frequency in this case) whereas 

the Timoshenko theory will give 


w(x. 


t) A cos (co t + c})) 
o ' o ^ 


+ B cos(w t + i/j) 
s 


(4-16) 


since it possesses two 

phase angles, and 

j 3 w V 
and — .) 


characteristic 
depend upon 


frequencies , 


and 


o s 

the initial conditions for 


(The 

w 


The idea that the Bernoulli-Euler theory (Equation 4-15) is included 
within the more exact, Timoshenko theory (Equation 4-16) is conveyed 
by the fact that the latter reduces to the former as the "amplitude" B 
tends to zero in Equation (4-16). Of course, the appropriateness of one 
theory or the other now depends upon the relative magnitudes of A and 
B. ° 

In fact, these "amplitudes" are actually functions of space, i.e.. 


A = A (x) 
o o 

B = B(x) 


and it was pointed out previously that the ratio is significant 

near the wavefront, where shear and rotary inertia are important. 
Furthermore, for "short" times i.e., on the order of 


T^ < t ^ T . 
ts ring 


(4-17) 
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the wavefront is still relatively near the plane of loading (x = 0) 

and short wavelength, high-frequency waves are still present near the 

origin to give difficulty for Bernoulli-Euler theory. Thus, it is 

clear that the deflection at the origin, namely w(0, t) , cannot coincide 

identically for the two theories, at least beginning from time t = 0 

and until time t = . 

ts 

At the later times, i.e., t >> T , the thickness- shear mode has 

u S 

undergone many vibration cycles and indeed the ratio ^^7 then 

approach zero at the origin (x = 0). In this case, the more exact 
Timoshenko theory may degenerate into Bernoulli-Euler theory and the 
displacement time histories coincide at v = 0, (i.e., at least well 
away from the wavefront) . 

This intuitive discussion can be put on a more rigorous basis by 
using the Fourier Transform (k, x) and Laplace Transforms, (s, t) in 
conjunction with the appropriate shell theories. Approximate results 
(i.e., asymptotic expansions) can be obtained by expanding the transforms 
as the Laplace parameter s i.e., time t 0, etc. Indeed, readers 

thoroughly familiar with elastic wave propagation may be able to provide 
references where such an approach has been carried out for the loading 


p(x, 0, t) = I 5(x) 6(t) 
o 

and geometry considered here. However, it was beyond the scope of the 
present contract to generate or extensively search for such analytical 
solutions. 


Additional Experimental Results 
for the Thick-Walled Shell 

The lack of quantitative agreement between theory and experiment 
for the thick-walled shell was somewhat disturbing, especially in view 
of the following items: 

• Previous holographic measurement of wave propagation (in 
beams and plates, Refs. 7, 8 ) had been successfully 
correlated with computer solutions . 

• The results presented herein (Section 3.0) for the thin- 
walled shell, agree reasonably well with calculations 
(e.g., see Figure 26). 

• The computer program (Ref. 17 ) has been amply verified by 
testing it against known, previous solutions, and finally 

• Examples are available (e.g.. Ref. 27 ) which show that (at 
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least for specific problems ) the solutions using Timoshenko 
shell theory approach those of Bemoulli-Euler theory 
(at sufficiently long wavelengths). 

In an effort to clarify and to better understand the thick shell 
results, additional interf erograms were analyzed to provide more 
deflection data. The experimental results, showing the radial deflection, 
w, vs. axial position, x, (for several discrete times) are shown in 
Figure 46. This family of curves illustrates typical wave propagation 
behavior, with the wavefront advancing (to the right) along the shell. 

This high-frequency, short-wavelength region occurs at the front of the 
wave, with the longer waves propagating more slowly and lagging behind. 

A feature of some interest is the behavior at the origin (x = 0) , as a 
function of time, t. The central def lection(s) w(0, t^) from Figure 46 
have been cross-plotted to give the time-history, w(0, t) , shown in 
Figure 47, The relatively smooth, continuous nature of the (albeit 
limited) experimental data in Figure 49 give strong support to the 
"consistency" of the present tests. 

Also shown in Figure 47 is the calculated half-period, T^/2 = 75 ysec, 

of the axi-symmetric ring mode of the shell. From the discussion of 
length scales and "characteristic times" presented previously, one might 
expect Bernoulli-Euler shell theory to predict a deflection time-history 
w(0, t) involving a half period of T^/2* On the other hand, Timoshenko 

theory will involve half-periods of T /2 and T^ /2, and combinations 

.-X. r o' ts * 

thereof. 

From an experimental standpoint, there are many more "characteristic 

times" in the physical problem than are indicated by Equations (4-12) 

and (4-13). For example, the input pulse has a duration T = 10 ysec, 

which may have a significant effect on the shell behavior. In addition, 

the finite width of the load, -e < x < e, of about 1/4 of an inch, 

introduces another characteristic length (for example, a flexural 

wavelength) which has an associated frequency, etc. Thus, it is not 

surprising to find an experimental "half-period" of approximately 50 ysec 

vs. an (idealized) result of T /2 = 75 ysec. 

o 

A thorough (analytical) study of this problem, including the effects 
of a finite width of the loading (in space) and a finite pulse duration 
(in time) would be of considerable interest, but such an analysis was 
beyond the scope of the present study. 
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Radial displacement, w, (inches) 



Figure 46: Radial deflection, w, vs. axial position, x. (Thick shell, 

no cut-outs.) 
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Central deflection, w(0, t) (inches) 






Results from Elasticity Theory 


Consider the following argument. The pulse duration, T = 10 psec, 
is a characteristic time, which can be used to calculate a structural 
wavelength, A . In other words. 


A = 2c T = 2(-l in/ysec) (10 ysec) 
s s 

A =2 inches, based upon the shear wave speed, c 
s s 

Having this value of A , and knowing the radius a =5-0 inches , 
the thickness h = .25 inches, we have 


h 

a 


■ 25 
5 


.05 


h 

A 


■ 25 
2 


.125 


With these values, we enter Ref. 28 and use the tabulated results for 
the n = 0 (axi- symmetric) mode. Page 20, of Ref. 28 then gives 


= .034 

for the corresponding frequency. But 


Q 


1 



.034 


7T C 

where o) = — r— = Fundamental thickness/shear frequency- 

5 O 


Thus 0 ) = -034 0 ) 

s 

and the corresponding period of vibration is 

X = ^ X L_ 

elast w 0 ) .034 

s 
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But 


2ir 

0) 

s 


T =4.1 iJsec, (see Section 3.0) and we then have 
s 


T - ^ = 4.1 ysec x 

elast .034 


= 120 ysec 


which results in a half-period of 


elas 


= 60 ysec 


which is very close to the experimental result (T/2 = 50 psec) shown 
in Figure 47. 

Of course, it is recognized that the duration of the input pulse, 
namely T = 10 ysec is actually 10 psec + (some fraction) which 
influences the structural wavelength A. Furthermore, the theory of 
Ref. 28 is calculated for the assumption of steady, traveling wave 
trains with wavelength A. Thus, it is considered to be good agreement 
when the experimental half period is 50 psec and the elasticity theory 
gives T/2 = 60 psec. 
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Test Results for the Thin-Walled Cylinder 


Wave Propagation Past a Circular Cut-Out - Experiments were 
conducted in which transverse waves were excited (and deflections 
recorded holographically) to study waves interacting with a circular 
cut-out. The thin-walled cylinder (h = .106 inches) was used in these 
tests, and the diameter of the cut-out was d = 2.0 inches. The impulsive 
load was applied at the longitudinal station x = 0, and the cut-out 
was centered at x = 6.0 (inches) to the right of the loading plane. 

Tests were conducted following the procedures outlined previously 
in Section 3.0. A typical inter ferogr am of the transverse wave inter- 
acting with the cut-out (at time t = 59 psec) is shown in Figure 48. 

The corresponding deflection data (w vs x) are presented in Figure 49 , 
for the centerline, i.e., 0 = 0. Non- symmetric behavior, caused by 
the interaction of an axi-symmetric wavefront with the cut-out, is 
evident in the interferogram. The high-frequency, short-wavelength, 
ripples (in the circumferential direction, 0, Figure 48) indicates 
the presence of circumferential harmonics (i.e., w cos n0) of high 
wave-number, n. 

Additional interf erograms , showing the wave passing the cut-out and 
the formation of a "shadow zone" to the right (beyond the cut-out) are 
presented in Figures 50 through 52. (For the earlier times, 
t ^ 59 ysec, see Appendix E.) The presence of short-wavelength, 
circumferential harmonics near the cut-out are evident in the interf erograms . 

Wave Propagation Past a Square Cut-Out - Similar tests were conducted 
on wave propagation past a square cut-out in the thin-walled shell. The 
square cut-out was 2.0 inches (on a side) and centered at x = -4 (inches). 
That is, the cut-out was centered 4 inches to the left of the loading 
plane, x = 0. Several interf erograms (which show the wave approaching 
and engulfing the cut-out) are presented in Figures 53 through 58. 

A typical plot of the radial deflection for this problem is given 

in Figure 59, t = 60 ysec. Note the large gradients (^) near the edge 

of the cut-out; see Figures 54 and 55. Non-axisymmetric behavior 
is apparent in the interf erograms , as indicated by the circumferential 
ripples in the dark interference fringes. The "shadow zone" caused by 
the cut-out is quite evident as well, c.f. Figure 58. Although it 
would be desirable to reduce these highly "two-dimensional" results 
to give quantitative data (e.g., in both x and 0) such a reduction 
was not possible during the present study. 

Additional results (showing the behavior at other times, t^) are 

given in Appendix E. When computer programs become available which can 
(economically) handle such non-symmetric. two-dimensional behavior, it is 
hoped that they will be compared with the holographic results given 
herein. 
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Radial Deflection, w. Inches 
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Figure 49: Radial Deflection, w, as a Function of Position along the Shell, 

with a Circular Cut-Out; see Figure 48). 


(Wave Interacting 


00 

00 



Figure 50: Transverse Wave Interaction with a Circular Cut-Out. 

Time t = 70 ysec 





I 



Figure 51 : Transverse Wave Interaction with a Circular Cut-Out. 

Time t = 80 ysec 
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Figure 53 


Wave propagation past a square cut-out 
h = .106 wall thickness.) 


Time t = 39 ysec. (Thin shell 
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Figure 54s Wave propagation past a square cut-out. Time t = 50 ysec. (Thin shell, h = .106 
wall thickness.) 
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Figure 56: Wave propagation past a square cut-out. Time t = 71 ysec. (Thin shell, h 

wall thickness.) 









Figure 57 : 


Wave propagation past a square cut-out. 
wall thickness.) 
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Figure 58: Wave propagation past a square cut-out. Time t = 92 ysec. (Thin shell, h =* 

wall thickness.) 
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Figure 59; Radial deflection^ vs. axial position, 
c.f. Figure 55.) 


h = .106 inch 
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Wave Interaction with a Stiffening Ring - The final test series 
(using the thin-walled shell) was designed to record the interaction of 
an axisynunetric wave with an axisymmetric ring stiffener* The ring 
stiffener had a rectangular cross-section with a width of .5 inches and 
a depth of 1.0 inches. It was centered approximately 4.0 inches to the 
left of the loading plane (x = 0). 

Typical interfere grams for this problem are shown in Figures 60 through 
63. The appearance of non-symmetric modes (i.e. cosn0) when the wave hits 
the stiffener are evident in Figure 62. Deflection data (w vs. x) for one 
interferogram are shown in Figure 64. Additional interfere grams are pre- 
sented in Appendix E. 
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Figure 60 : Wave Interacting with the (Internal) Stiffening Ring 

(on the left, at x = -4 inches) and the Cut-Out (on 
the right). Time t = 30 ysec 
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Figure 62 : Wave Interacting with the (Internal) Stiffening Ring 

(on the left, at x = -4 inches) and the Cut-Out (on 
the right) . Time t = 49 usee 
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Figure 63: Wave interacting with the (internal) stiffening ring (on the left, 

at X = -4 inches) and the cut-out (on the right). Time t = 60 ysec. 
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Figure 64: Radial Deflection, w, as a Function of Position along the shell. (Wave 

Propagating to the left, past the Stiffening Ring at x = -4.) 



5.0 WAVE. PROPAGATION IN CONICAL SHELLS 


Introduction and Background 

In addition to the tests on cylindrical shells, it was desired 
to obtain data for wave propagation in conical shells as well. Two 
conical specimens were tested, and the results of the tests are given 
in the following paragraphs. 


Tests on the Cone-Cylinder Specimen 

A thin-walled cone- cylinder specimen was obtained for use in 
these tests. The dimensions of the specimen are given in Figure 65, 
and a photograph of the actual model is given in Figure 66 . 

In order to apply the loading pulse (electromagnetically) the 
loading wire was located in the larger end of the specimen. Initial 
tests were conducted with the load applied at x = 1.0, i.e., approximately 
1 inch from the large end of the cone. 

A typical interf erogram for this test arrangement is given in 
Figure 67. The corresponding deflection data are given in Figure 68. 

As the waves propagate longitudinally along the cone, the appearance 
of non-symmetric modes (i.e., cos n 0) becomes evident. The development 
of the non-symmetric response is evident in Figures 69 through 72. 
Additional interferograms are given in Appendix F. 


Tests on the Large Conical Shell 

In order to investigate waves propagating in a larger shell, 
another conical specimen was built. A-sketch showing the shell 
geometry is given in Figure 73. A photograph of the shell itself 
is given in Figure 74. 

The test procedure and experimental set-up were similar to those 
discussed previously. The shell has a thickness h = .094 inches, 
and it was loaded with the repelling-wire device discussed previously. 
A typical hologram for this test series is given in Figure 75 . The 
corresponding deflection w (as a function of the length, s, along 
the cone) is given in Figure 76. 

Additional interferograms are given in Figures 77 through 79, 
which again demonstrate non-symmetric behavior (i.e., cos n 6 modes). 
Finally, additional interferograms are presented in Appendix F. 
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NOTE: 

1. Circular Section 

2. Aluminum 


Figure 65: Dimensions for the thin-walled cone cylinder. 




Figure 66: Photograph of the thin-walled cone-cylinder (diameter 

of large end is 2.6 inches). 
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-walled cone-cylinder. Time t = 10 psec. 




Radial Deflection, w, (inches) 



Figure 68: Radial deflection, w, vs. axial position, x. (Load centered at 

X = 0.) Cone cylinder combination, time t = 10 -sec. 


(inches) 
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Figure 69: Wave propagating in a thin-walled cone-cylinder. Time t = 48 ysec. 
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Figure 70 : Wave propagating in a thin-walled cone-cylinder. Time t = 51 ysec. 















Figure 74: 


Conical shell specimen (14° half-angle, 
conical f rustrun) . 
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Figure 75: Wave propagation in a conical shell. Time t = 29 ysec 







Transverse deflection, w, (inches) 






Figure 73 : Wave propagating in a large conical shell. Time t = 48 psec. 
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6.0 CONCLUDING REMARKS 


In an attempt to extract the major conclusions found in this study, 
the following points stand out: 

o Transverse wave propagation in shells can be accurately 
measured using pulsed laser holography. From an engi- 
neering standpoint, this means that holography can provide 
engineering data on the effects of cut-outs, discontinuities, 
etc. 

o The axi-symmetric response of a circular shell to an axi- 
symmetric loading pulse may subsequently become unstable 
and excite circumferential harmonics in the shell. This 
is a problem in dynamic stability, and may be worth 
investigating analytically. The presence of non-symmetric 
modes might be of practical importance since they will 
affect stress distributions, circumferential bending, etc. 

o When the input pulse is of sufficiently short duration, it 
is possible to excite responses that are beyond the scope 
of conventional shell theories. In this case, it is 
necessary to use "elasticity-type" solutions, as discussed 
in Section 4.0. The reader should note that elementary 
shell theory (Bemoulli-Euler) overpredicted the thick- 
shell deflections (as compared with the experimental data). 

Thus it seems that the elementary theory is overconservative 
(at least for this problem) and might result in a non- 
optimum design. 

o Holography can be used to measure waves propagating past 
cut-outs and/or stiffening rings, which exhibit pronounced 
non-symmetric responses. Dynamic stress concentration 
caused by cut-outs, discontinuities, bolt holes, etc., 
can be accurately evaluated using pulsed laser holography. 

Such stress concentration factors often have a major impact 
on the design of shell structures. 
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APPENDIX A 


THE EXPERIMENTAL LOADING TECHNIQUE 


Introduction 

In loading the shell (s), the objective was to produce a radial pressure 
that was uniform around the circumference, concentrated at a single axial 
station (x = 0) and of short duration in time (i.e., on the order of a few 
microseconds). The loading was to generate axis 3 nmnetric, transverse waves 
in the shell(s), which involve the radial displacement, w, and the longi- 
tudinal displacement, u. From the standpoint of analysis, the desired 
loading might be idealized as 


p(x,e,t) = P(t)6(x) 


(A-1) 


where d(x) is a Dirac delta function in space and P(t) represents the time- 
history of the loading pulse. In practice, P(t) can be approximated by 


P(t) 


A sin 


2 7Tt 
T 


0 ^ t <_ T 

(A-2) 


0 


T < t < 


where T is the time duration of the pulse. 

In launch vehicle applications, the loading function p(x,6,t) usually 
results from pyrotechnic shock loads (i.e., explosives). Thus, one of the 
first loading techniques investigated was a sprayed-on, light-sensitive 
explosive (Ref. 29). Such explosives have been used previously in structural 
applications, but we experienced two main difficulties with this technique: 

6 2 

(1) To cause detonation, high-intensity light (i.e., 7 x 10 lumens/in ) 
must be applied over all areas covered by the explosive. 

(2) A long turn-around time is required between tests, which involves 
spraying the explosive, drying it, etc. 

As a means of detonating the explosive simultaneously (around the cir- 
cumference of the shell) an exploding wire (Ref. 30) was considered. The 
possibility of loading the shell directly (by an exploding wire) led in turn 
to the development of the final loading device, which involves an electro- 
magnetic impulse. 
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The Repelling-Wire Loading Device 


Consider two parallel, current-carrying conductors, a distance d apart 
(see Fig. A,l). If one conductor carries a current Ij^(t) in one direction. 


i 


Figure A.l: Parallel, current-carrying conductors, 

snaced a distance d anart 


and the other conductor carries l 2 (t) in the other direction, then the force 
(per unit length) between the conductors is given by (Ref. 31) 


F(t) 


2-nd 


(A-3) 


where = 4ttx10 in the rationalized MKS system. F(t) is force of magnetic 
repulsion, since the interacting electromagnetic fields of the conductors 
cause them to fly apart (when and I 2 are in opposite directions). If one 

conductor is held fixed in space and the other is nestled against the 
structure to be loaded, it is possible to generate pressure- time histories 
of the desired magnitude and duration. 

To initiate waves in our circular shell structures, a flat aluminum 
conductor was formed in the shape of two concentric circles (see Fig. A. 2). 

The conductor did not form a complete circle, since it was necessary to 
bring the current in at one point on the circle and take it out at another. 
Thus, the circular conductor was made of one continuous piece of aluminum 
(1/4- inch wide) looped back upon itself. A rubber tube was installed 
between the conducting loops of aluminum; the tube and conductors were held 
in place by a circular plexiglas disc (see Fig. A. 3). The rubber tube was 
pressurized to force the loading wire into close contact with the shell. 

The generation of significant force pulses with a repelling wire device 
requires very high electrical currents, e.g., on the order of 50,000 amps. 
These currents are produced by charging a capacitor bank to several kilovolts , 
and then discharging the capacitors through the repelling wire circuit. The 
system can be adequately represented by a simple R-L-C circuit, as discussed 
in Ref. 32. (A circuit diagram is given in Figure A. 4). 
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Wire 

Spacing 


Figure A. 2: Concentric current-carrying loop(s) used to load 

the shell. 
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Figure A. 5(a): Current vs. time in repelling wire circuit 

(charging voltage 12 kv, sweep speed 5 ys/cm, 
scale 50 mv/div) . 
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Figure A. 5(b): Current vs. time in repelling wire circuit 

(charging voltage 20 kv, sweep speed 5 ys/cm, 
scale 50 mv/div) . 


Current in the Repelling Wire Circuit 
The response of the rep el ling-wire circuit is given by (Ref. 32) 


I(t) = — e sxnwt 
ojL 


(A-4) 


where 


0) 



4l2/ 


a 


2L 


is the circular frequency 
is the damping 


and is the initial charging voltage. 

Equation (A-4) refers to an underdamped circuit, and a typical oscillo- 
scope trace (current vs. time) for the experiment is given in Figure A. 5. 

The damping, R, was adjusted experimentally to yield the desired current- 
time history. From the trace of Figure A. 5, the circuit inductance L was 
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calculated, using the values C = 7 x 10 farads (of the capacitor bank) 
and R =• .53 ohms (measured with a micro-^ohmmeter) , and o) = 7 t/16x 10*’^ sec. 
Then the current, I(t) is given by 




-.0883t 

e 


sin(.2t) 


(A-5) 


where t is in microseconds . 


Calculation of the Force and Impulse 
Using Equations (A-3) and (A-4) , the force/length is given by 


F(t) = 


-7 2 

2x10 




2x7 f^oY -2at . 2 

e sin wt 


(A-6) 


The (impulse/unit length) that the wire exerts on the shell is obtained 
by integration: 


-X 


^ final 


F(t)dt 


(A-7) 


If is taken as 16 psec (the first half-period of the pulse) , then 


t = 16 Msec 


I 

o 


f -2at 
= J e 
o 


2 1 
sin tot dt = -'[.409] 

to 


(A-8) 


If the upper limit is taken as infinity, (which adds up all the ever-decreasing 
pulses) then 


t = “ 

= J* Sin^ojt dt = •^[•475] (A-9) 

o 


Using the experimental values for the charging voltage, V , and the wire 
spacing, d, we have 
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V = 19,000 volts 
o ’ 


I = 
o 


“2 

d = .45 in. = 1.14 x 10 meters 


/l.9 X 10^ \ 2 X 10 / 1 \ /.409 \ 

' / \2 X 10^ / \1.14 X 10~V 


^o- {4’.2o} "lO 


-'2 Newton- sec 
meter 


(A-10) 


where the two values given depend upon the upper limit taken for the integral 
(see A-8 and A-9). 

In English units, the specific impulse is 


J2.06( 
(2.39 f 


X 10 


-4 lb-sec 
in 


(A-11) 


using equation (A-10) and appropriate conversion factors. 


Measurement of the Impulse Using A Ballistic Pendulum 


An attempt was made to measure the specific impulse (I ) directly, using 
a ballistic pendulum. For this purpose, the repelling wire was installed in 
a calibration ring (see Figure A. 6). The ring was of the same diameter, 
material, and thickness as the shell to be tested. A small square hole was 
cut through the wall of the calibration ring to accommodate a small ballistic 
pendulum. The pendulum mass was a small square piece of aluminum, which was 
suspended on long strings and then nestled against the repelling wire. 

The capacitor bank was charged to a known voltage, V , and then the 
switch was closed to ”fire" the circuit. The deflection of the pendulxmi was 
measured, and the impulse values (per unit length) were calculated. A plot 
of Iq (specific impulse) vs Vq (charging voltage) is given in Figure A. 7, 
for a wire spacing of d = .45 inches. 

At 19KV, the measured impulse was 


(I ) = 1.34 X 10 

o meas 


-4 lb-sec 


in 


(A-12) 


vs a calculated value of 
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Ji 





calc 


2.39 X 10 


lb-sec 

in 


from equation (A-11) . 

As things turned out, the experimental deflection data compare more 
favorably with the analysis when the larger value (2.39 x 10 ) is used in 

the computations. This led to the conclusion that the ballistic pendulum 
technique does not correctly measure the specific impulse, I , imparted to 
the shell. Apparently, the short pendulum mass does not respond the same 
as the continuous shell. This anomalous behavior is attributed to eddy 
currents thought to be present in the continuous shell but largely 
lacking in the pendulum mass. 


Description of the Experimental Pulse in 
a Form Convenient for Analysis 

As presented in Section 3.0, it is convenient to represent the 
applied loading, p(x, 0, t) , as 


p(x, 6, t) = A sin^ ^ (A-13) 

within the domain 

-e < X < e 

0 < e < 2tt (A-14) 

0 < t < T 

and zero everywhere else. 

The x-domain (with e = 1/8 of an inch) 


- 1/8 < X < 1/8 

in Equation (A-13) is based upon the fact that the repelling wire (which 
loaded the shell) was a flat, aluminum strip, 1/4 of an inch wide. No 
attempt was made to measure the variation of the applied pressure loading 
in the x-direction. In fact, a more realistic approximation for the 
x-dependence might be 


within 


2 

p(x, 0, t) cos 4t7X 
-1/8 < X < 1/8 


(A-15) 
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Similarly, the exact variation of the pressure in the circumferential 
direction was not measured. The loading wire which applied the pressure 
was not a complete circle around the circumference. It was necessary 
to bend the wire (to accommodate electrical leads, etc.) in such a way 
that a gap existed in the 0-directlon. (For example, see Figure A. 2.) 

A more appropriate description of the loading (in the circumferential 
direction) might thus be given by 


p(x, 0, t) A^ cos n0 (A-16) 

n=0 

where the axi-symmetric term (A ) is much larger than the circumferential 
harmonics (A ) : 


A >>A n=l,2,3, .... 

o n 

Finally, the variation in time is more exactly given by Equation 

(A-6) 


F(t) = e sin^ cut 
where a = .0883 x 10^ sec ^ 

and 03 = . 2 X 10^ sec ^ 

Thus it is apparent that the experimental loading approximates 
(but by no means duplicates ) the form given by Equation (A-13) and 
used in the analysis. 
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APPENDIX B 


WAVES IN THE THIN-WALLED CYLINDER 
(ADDITIONAL INTERFEROGRAMS) 


This Appendix is meant to supplement Section 3.0 of the report, and 
it contains results for several times, t^. 
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Figure B.2: 


Wave propagation in the thin-walled cylindrical shell. Time t 


99 usee. 
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Figure B.4; Wave propagation in the thin-walled cylindrical shell. Time t = 150 psec. 
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Figure B,5: Wave propagation in the thin-walled cylindrical shell. Time t = 171 ysec. 
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Figure b. 6^ Wave propagation in the thin-walled cylindrical shell. Time t = 197 ysec. 








APPENDIX C 


WAVES REFLECTING FROM THE ENDS 
OF THE SHELL 

(ADDITIONAL INTERFEROGRAMS) 


The interfero grains presented in this appendix show transverse 
waves reflecting from the ends of the thin cylindrical shell, and they 
illustrate the structural behavior at various times, t^. These figures 

are intended to supplement Figures 20, 22, and 24 in the main, body of 
the text. Section 3.0. 
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Figurec,2* Wave Propagating Toward the "Clamped" End (at the Left), 
Time t = 30 psec. 



o 



Figurec.3: Wave Propagating Toward the "Clamped" End (at the Left). 

Time t = 60 ysec. 
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FigureC.4: Wave Propagating toward the "Clamped" End (at the Left). 

Time t = 75 ysec. 
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Figure c. 6: 


Wave Propagating toward the "Clamped” End (at the Left) 
Time t = 105 ysec. 
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Figure C.& 


Wave Propagating toward the Free End (at the Right) 
Time t = 75 usee. 
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Figure c. 10: Wave Propagating toward the Simply-Supported Boundary 
(Vee-Groove) at the Left. Time t = 11 ysec. 
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Figure c, 12 : Propagating Toward the Simply-Supported Boundary 

(Vee-Groove) at the Left. Time t = 30 ysec. 






I 
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Figure C. 14- Wave Propagating Toward the Simply-Supported Boundary 
(Vee-Groove) at the Left. Time t = 48 psec. 




159 



Figure c. 16: Wave Propagating Toward the Simply-Supported Boundary 
(Vee-Groove) at the Left. Time t = 71 ysec. 


APPENDIX D 


This 
Section 4 


WAVE PROPAGATION IN THE THICK-WALLED SHELL 
(ADDITIONAL INTERFEROGEAMS) 


Appendix is meant to supplanent the results given in 
0, "Tests on a Thick-Walled Cylinder." 
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Figure d,i: Wave propagation in the thick-walled shell, (h = .25 inches). 

Time t = 20 ysec. 


Figure d, 2: Wave propagation in the thick-walled shell. (h = .25 inches) 

Time t = 30 usee. 
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Figure d. 3* Wave propagation in the thick-walled shell. (h = .25 inches). 
Time t = 47 psec. 
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Figure D.4 i Wave propagation in the thick-walled shell. Time t = 51 ysec. 
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Figure D. 5 : Wave propagation in the thick-walled shell. Time t = 53 ysec. 




Figure D.6: Wave propagation in the thick-walled shell. Time t - 54 ysec. 
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Figure d. 9* Wave propagation in the thick-walled shell. Time t = 101 psec. 
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Figure D.ll: 


Wave interacting with cut-outs. Thick-walled shell, 
h = .25 Inches. Time t = 40 vsec. 
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APPENDIX E 


WAVES PAST CUT-OUTS AND STIFFENERS 
(ADDITIONAL INTERFEROGRAMS) 


This Appendix is to supplement Section 4.0, in particular, “Waves 
in the Thin-Walled Cylinder.** 
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Figure e.3 


I 


Transverse Wave Interaction with a Circular Cut-Out, 
Time t = 30 ysec 
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Figure E.7- Wave propagation past a square cut-out. Time t 
h “ .106 wall thickness.) 


11 ysec. (Thin shell, 
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Figure E^IO: Wave Interacting with the (Internal) Stiffening Ring 
(on the left, at x = -4 inches) and the Cut-Out (on 
the right) . Time t = 10 psec 






Figure E.ll*- Wave Interacting with the (Internal) Stiffening Ring (on 

the left, at X = -4 inches) and the Cut-Out (on the right). 
Time t = 21 ysec 
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Figure E. 13: Wave Interacting with the (Internal) Stiffening Ring (on 
the left, at x =-4 inches) and the Cut-Out (on the right). 
Time t = 80 ysec 
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APPENDIX F 

WAVES IN CONICAL SHELLS 
(ADDITIONAL INTERFEROGRAMS) 


This Appendix is designed to supplement Section 5.0 of the main body 
of the text. 
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Figure F.2; Wave propagating in a thin-walled cone-cylinder. Time t = 12 ysec. 
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Figure 5 Wave propagating in a thin-walled cone-cylinder. Time t = 28 psec. 
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Figure F.8: Wave propagating in thin-walled cone-cylinder . Time t = 10 psec. 
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Figure F.IO: Wave propaga 


in thin-walled cone-cylinder. Time t = 30 ysec. 
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Figure F.13: Wave Propagating in a large conical shell. Time t = 20 ysec. 
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Figure F.15: Wave propagating in a large conical shell. Time t = 81 ysec 










